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Abstract: In the presence of compact dimensions massive solutions of General Relativity 
may take one of several forms including the black-hole and the black-string, the simplest 
relevant background being M 3+1 x S 1 . It is shown how Morse theory places constraints 
on the qualitative features of the phase diagram, and a minimalistic diagram is suggested 
which describes a first order transition whose only stable phases are the uniform string and 
the black-hole. The diagram calls for a topology changing "merger" transition in which 
the black-hole evolves continuously into an unstable black-string phase. As evidence a 
local model for the transition is presented in which the cone over S 2 x S 2 plays a central 
role. Horizon cusps do not appear as precursors to black hole merger. A generalization to 
higher dimensions finds that whereas the cone has a tachyon function for d = 5, its stability 
depends interestingly on the dimension - it is unstable for d < 10, and stable for d > 10. 
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To my wife, 
Dorit. 

1. Introduction 

In a gravitational background with compact dimensions several phases of black object so- 
lutions to General Relativity may exist, some of them having distinct horizon topologies, 
depending on the relative size of the horizon (say the Schwarzschild radius) and relevant 
length scales in the compact dimensions. As the mass of the object is varied phase tran- 
sitions are expected to occur between the various phases. For concreteness, consider the 
simplest case -3+1 extended dimensions and a fifth compact dimension of radius L, which 
is denoted here by the z coordinate. Black objects with Schwarzschild radius much smaller 
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than L are expected to closely resemble a 5d Schwarzschild black hole (BH) with a nearly 
round S 3 horizon topology, since by the equivalence principle local measurements at the 
vicinity of the BH should not detect the large finite size of the z axis. Very massive black 
objects on the other hand, will have a horizon size much larger than L and will be extended 
over the z axis wrapping it completely. Such objects will have S 2 x S 1 horizon topology 
and will be referred to as "black strings" . One type of black string is the 4d Schwarzschild 
solution with the z coordinate added in as a spectator and it will be called here a "uni- 
form string", while another type of solution may be z dependent, namely, a "non-uniform 
string" . 

An attempt to build a consistent phase diagram for this system will lead us to expect 
a novel topology change within classical General Relativity. The transition involves a one 
parameter family of Ricci-flat metrics which passes through the cone over S 2 x S 2 (locally). 
This cone may be deformed in one of two ways much like the conifold transition leading 
to different topologies. Each deformation keeps one incontractible sphere, and "fills-in" 
the other, so the transition may be described as an S 2 shrinking to zero size and then a 
different S 2 growing up. In addition to the similarities with the conifold the differences 
should be pointed out as well - this transition happens in any dimension higher or equal 
to 5, and there are no Killing spinors. 

Gregory and Laflamme (GL) |l|, |2| found that the uniform black string solution develops 
a perturbative instability at a certain critical radius, which they interpreted as a transition 
to the black hole. Horowitz and Maeda || recently surprised many by arguing that the 
endpoint of the decay of the string could not be the black hole, definitely not in finite 
affine parameter, basically because of the "no tear" property of horizons. Later Gubser 
H showed that the transition is actually first order by studying local properties of the 
Gregory-Laflamme point. 

These papers studied the transition in one direction, namely starting with the black 
string and reducing its mass, which will be denoted here by "string — > BH" (not implying 
that the end point is necessarily a BH). One may consider also the opposite direction - 
starting with a 5d BH and increasing its mass - this direction will be denoted by "BH 
— > string". Regarding this latter direction, Susskind |J claimed (based on Matrix M 
theory, analogy and intuition) that the transition from a black hole to a non-uniform 
string (different from the GL point) would be analogous to a Gross- Witten transition in 2d 
gauge theory, a third order transition (a transition so smooth that only the third derivative 
of the free energy jumps). 

The original goal of this research was be to determine the phase diagram for this 
transition, and especially take the "BH — * string" route and identify the critical point 
where the black hole develops a perturbative instability, expecting it to be related to the 
self intersection of the horizon - the merger of the north and south poles across the compact 
dimension. 

The motivations were 

• A case study of a rich GR phase diagram. 

• A system with a horizon merger, without time dependence, which could be a model 
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(both analytically and numerically) for the heavily studied problem of black hole 
collision. 

In particular, numerical simulations in 4d show (see the review [[?J and references 
therein) that when two BH's approach each other a cusp forms on the horizon of 
each BH facing the other one. Would cusps be present here? 

• A test for the controversy raised by ||. 

One would like to determine whether there exists a third non-uniform phase as 
claimed by ||, or features of Gross- Witten transition as claimed in §. 

• A challenge for a quantum interpretation of a topology change in the horizon. 

In string theory or any other theory of quantum gravity a black hole is a high entropy 
quantum system (for example, the black- hole D-brane model). The process of black 
hole merger should translate to increased coupling between the respective quantum 
systems. It would be interesting to understand the irreversible merger point in the 
quantum theory together with the continuous evolution of the entropy as the two 
BH's approach each other, but it will not be discussed in this work. 

• Wide relevance for higher dimensional branes and extra dimensions. 

In string theory branes and particles in low dimension often can be interpreted as 
some fundamental branes in lOd or lid, and then the issue arises whether they should 
be considered localized or smeared over the compact directions, which is exactly the 
transition studied here. This system shows up in many other contexts in string theory 
as well. 

Other contributions include the realization of the connection between the GL instability 
and negative specific heat P, 01, a discussion of the charged non-uniform black string 



|1C], an attempt to use the Petrov classification to derive the metric of the non-uniform 
string pH , a numerical study of a closely related problem and other recent papers 
|T|, [15], 0, 0, 0,0, §0|. 



We start in section |2| by setting up the system and our notation and reviewing a number 
of central contributions which we summarize by an incomplete phase diagram (figure |||), 
looking like an archeological finding with missing parts. 

In section |3| Morse theory is used (after a short review) to demonstrate important 
constraints on the phase diagram. I argue why these constraints require a topology change 
transition and disfavor a stable non-uniform phase. Finally I present the simplest (in my 
opinion), consistent phase diagram (figure [R]). 

In order to find evidence for a qualitative picture, more qualitative tools (beyond the 
Morse index) are sought. In section |] the topology of the various solutions is analyzed, and 
it is found that the suggested phase diagram requires a topology change - a continuous 
path of solutions connecting different topologies. That leads me to seek a local model of 
the topology change within Ricci-flat metrics. Moreover, I conjecture that the point of 
topology change is also a phase transition point. In the rest of the section, the topological 
requirements of the local transition are elaborated. 
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The expected Ricci-flat solutions are sought for in sections || ^ and a final summary 
and discussion is presented in section [?]. 

Note added: This paper is the published version of [ 51 1 . Only minor phrasing changes 
were made to improve the text but it was not updated. For updates see the review [52]. 



2. An incomplete phase diagram 

Consider a static (no angular momentum) black object in an M 3+1 x S 1 gravitational 
background, namely extended 3 + 1 space-time with a periodic fifth dimension which will 
be denoted by the z coordinate. The system is characterized by 3 dimensionful constants: 
L the size of the extra dimension z ~ z + L , M the (4d) mass of the system measured 
at infinity of M 3 , and G5 the 5d Newton constant. From these a single dimensionless 
parameter can be constructed 

p = (G 5 M)/L 2 , (2.1) 

while the 4d effective Newton constant is given by G4 = G5/L. 

The isometries of these solutions are SO(3) x U(l), where the SO(3) comes from the 
spherical symmetry in M 3 and the U(l) comes from time independence. 1 The most general 
metric with these isometries is 

ds 2 = -e 2A dt 2 + ds\ r>z) + e 2C dfl 2 , (2.2) 

which is a general metric in the (r, z) plane together with two functions on the plane 
A = A(r,z), C = C(r,z). The horizon is a line determined by e 2A = and d£l 2 = 
d9 2 + sin 2 6 dip 2 . In 4d one can choose a coordinate system that is especially adapted to 
cylindrical symmetry |2l}| , but this works best only in 4d. 

At least two phases of solutions can be distinguished - a black string with an S 2 x S 1 
horizon topology and a 5d black hole with an S 3 horizon. The uniform black string (figure 
[[]) is given by the 4d Schwarzschild metric with z added as a spectator coordinate 

ds 2 = -(1 - r 4 /r) dt 2 + (1 - r 4 /r)~ 1 dr 2 + dz 2 + r 2 dCl 2 (2.3) 

where the Schwarzschild radius is r 4 = 2 G4 M. 

The small 5d black hole (see figure |)(r 5 = ((8/3vr) G 5 M) 1/2 < L) can be approxi- 
mated by a combination of two solutions. Denoting the distance from the black hole by p, 
for p <C 7"5 it can be approximated by the 5d Schwarzschild BH 

ds 2 = -(1 - r 5 2 /p 2 ) dt 2 + (1 - r 5 2 /p 2 y l dp 2 + p 2 dn 2 s3 (2.4) 

where dQ, 2 sS = dx 2 + sin 2 ^^ 2 - For p 3> the Newtonian approximation is valid, and 
the potential is proportional to the Green function 

V(w) = — [coth(u>/2) + coth(u>/2)]. (2.5) 

1 In the Lorentzian solutions it is the non-compact version of U(l). 
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Figure 1: The uniform black-string, r is the radial coordinate in 



where w := 2 it (r + iz)/L. Indeed close to the source point w — > 0, the potential has the 
expected 5d behavior V ~ (l/4r)[(2/iy) + (2/w)] = l/(r 2 + z 2 ), while for Re(u>) — > oo, the 
4d behavior is restored coth(io) = 1 + 2exp(— 2w), V ~ (l/2r) + 0(exp(— r). Presumably 
a full solution can be built perturbatively in r§/L from these two approximations. Note 
that following the equipotential surfaces of (^^) one encounters already a primitive version 
of a topology change - close to the source the surface has an S 3 topology, then there is one 
singular surface with a conic singularity (a cone over S 2 ) after which the topology of the 
surfaces changes to S 2 x S 1 . 

Comparing the areas of the two solutions 
one finds that for the string the area is 
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while for the small BH 
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(2.6) 



(2.7) 
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Hence, for small \x the black-hole is preferred, 
while for larger \x it is not clear if the black- 
hole phase exists, but even if it does than if 
its area would grow as in ( |2.7D it would be 
dominated by the black-string. Therefore a 
phase transition is expected from the outset. 

Let us start by assembling the existing 
knowledge into an incomplete phase diagram 
(figure |3|) . The vertical axis is the parameter 
fj,, while the horizontal axis is an order pa- 
rameter 2 of non- uniformity A. For small A it is defined by A : — (1 — ^"min/ ^) where T m i n , v 
are the minimum and average respectively of r = r(z), so that A = exactly when the 
solution is a uniform string. For large A no precise definition is offered, but consider it to 
continue to be a measure of non-uniformity even in the black-hole phase. 



Figure 2: The 5d black- hole with Newtonian 
equipotential surfaces. 



2 Even though normally a phase diagram does not include order parameters. 
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Figure 3: Collecting present data into a phase diagram, /i is the dimensionless parameter of 
the problem, and A is an order parameter, a measure of non-uniformity along the z axis. Stable 
(unstable) phases are denoted by a solid (dashed) lines (double dashed lines denote phases with 2 
negative modes). GL denotes the Gregory-Laflamme critical point at \xql and GL' is a copy at 
Mgl/2. 



The horizontal axis is inhabited by the uniform string phase, which changes behavior 
around the numerically determined Gregory-Laflamme critical point ^ugl — -070 corre- 
sponding to (L/r^cL — 7.2 [|[. We denote the stable (unstable) phase for /i > fiQL 
(/i < /xgl) by a solid (dashed) line. The black-hole phase, drawn to the left, is known to 
exist for ji <C 1 but it is not clear how is evolves for fj, ~ 1. 

The non- uniform phase (figure ||) , argued by Horowitz and Maeda, is depicted for some 
intermediate values of \i where the decaying uniform string could end up. However, we 
must pause and be critical here. Since the predicted solution is not available yet, neither 
analytically nor numerically (|13, 111]), 3 we stay on the side of caution and mark it by a 



3 Note that |H| came out about a year ago. 
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question mark. 



Figure 4: The Gregory-Laflamme unstable mode, and a non-uniform black-string 



At the GL critical point the marginally tachyonic mode leads the way to a new phase. 
A priori this phase could be either stable or unstable (the two possibilities are shown in 
figure H). In a beautiful work, initiated by the hope of finding the stable non-uniform 
phase, Gubser concluded that the emanating phase is actually unstable ||. This was done 
by following the solution in a perturbative series near the GL point (to third order in 
some quantities!). Here again we must be critical, noting that || contains a so called 
"unwanted scheme dependence", namely the results depend on an unphysical parameter 
(this dependence looks to me like a missing boundary condition; the local analysis around 
the GL point is discussed further in the appendix). However, since the results are in 
some sense only weakly uncertain and since the bottom line agrees so far with unpublished 
numerical simulations [13, 14], we have some confidence in it. 

It should be mentioned that the GL critical point has infinitely many copies at fi = 
^Gh/n. These are simply a consequence of the unstable phase of equally spaced n black 
holes. Since all the phases connected to these copies are duplicates of the basic phases they 
will be altogether ignored in the rest of the paper. For example, at fi = a string 

with one unstable mode turns into a string with 2 unstable modes, and a new phase with 
2 negative modes emanates from the point. 

The phase transition, like any other one, can be approached from two different direc- 
tions, showing possibly different behavior. The two directions will be denoted by "string — > 
BH" and "BH — > string" but not implying the absence of other phases. So far we discussed 
only the string —* BH, but this work will concentrate on the opposite one. 
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Regarding the latter, BH — * string direc- 
tion (figure |5|), one should mention the 
ideas of Susskind | on a relation with 
the Gross- Witten transition in gauge the- 
ory [||], stemming from a matrix theory 
description. In the Gross- Witten tran- 
sition in 2d large N gauge theories one 
studies the distribution of eigenvalues of 
a unitary matrix, exp(i#), the matrix 
being the Wilson loop over an elemen- 
tary plaquette. The behaviour of the 
eigenvalues is given precisely by a me- 
chanical model of N point particles mov- 
ing on a vertical ring (the bottom of the 
ring is 6 = 0) in the presence of a uniform gravitational acceleration g, a repulsive potential 
logarithmic at small separations due to the Van-Der- Monde determinant and a temperature 
T (see figure ^ . It is found that one can take a large N limit in such a way that the nor- 
malized eigenvalue density distribution approaches a limit p{6) and a single dimensionless 
parameter remains, which is essentially T. Moreover, whereas the probability density at 
the top p(ir) never vanishes for finite N, it does vanish in the large N limit for small enough 
T. On the other hand, for very large T, the distribution becomes uniform p{9) = 1/(2 tt). 
The phase transition occurs at the critical temperature above which p(ir) > (see figure 
|6|), and it is found to be third order (so smooth that only the third derivative of the free 
energy jumps). 



Figure 5: Approaching merger. Do cusps appear? 



■JO 



Figure 6: The behavior of the gas of eigenvalues in the Gross- Witten phase transition. The three 
different density curves are the low temperature "blob" phase analogous to the black-hole, the 
critical curve, and the high temperature "winding" phase analogous to the black string. 

According to H, the distribution of eigenvalues along the 9 circle is analogous with 
the mass distribution of the black hole along the z circle, the low T "blob" of eigenvalues is 
analogous with the black hole phase, and the high T "winding" phase is analogous with the 
(non-uniform) string. ||] goes on to predict the existence of a Gross- Witten like transition 
in 4d large N gauge theory, based on the behavior of black holes in a background with a 
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compact 3-torus. 



3. Phase diagram - novel constraints 
3.1 Euclidean version and Free energy 

For Concreteness and convenience, we choose to use the Euclidean signature for the geome- 
tries, and to do thermodynamics with the free energy, namely in the canonical ensemble. 
Since the geometry is time independent we are free to use the Euclidean signature over 
the Lorentzian. In the Euclidean version "time" is periodic and the circle fiber degener- 
ates over the horizon, so the geometry cannot be continued inside the horizon. Another 
difference is that Euclidean black hole solutions have a negative mode not present in the 
Lorentzian version |^2| (which is actually the same as the GL mode, and its eigenvalue 
gives ugl)- Hence we make the convention that the number of negative modes indicated 
in the phase diagrams is 1 less than that number for the Euclidean version, to conform 
with the Lorentzian picture. From here on we avoid the Lorentzian signature and any time 
dependence until the discussion section. 

The canonical ensemble offers a direct expression for the free energy F = F{(3) through 
the action / 

-PF = I=—^—[ R+——f (K-K ) (3.1) 



16itGJm 8itG jgj^i 
where the first integral is the bulk contribution of the Ricci scalar and the second is a 
boundary contribution where K is the trace of the second fundamental form on the bound- 
ary, and Kq is the same quantity for a reference geometry which is flat K 3 x S 1 x S 1 in 
our case. When passing to the canonical ensemble (3 replaces M as a fundamental variable 
and so fi is redefined by 

"=871 < 3 ' 2 > 
where the constant factor was determined so that the new and old definitions would coincide 
for the uniform string. No confusion should arise due to the two definitions of [i Q2.1| , |3.2| ) 
since they coincide for the uniform string which was the only place we used \i so far, and 
from now only fj, = /x(/3) ( fT^ ) will be used. 

In our case, the free energy should be thought to be an even function of A. The trans- 
formation A — > — A is a symmetry of the action, a remnant of the symmetry of translations 
in the z direction, which was fixed by the symmetry breaking due to the GL mode. 

3.2 Morse theory 

Since our goal of understanding the structure of the phase diagram is qualitative rather 
than quantitative 4 we should seek appropriate qualitative tools. As we follow a curve of 
solutions, namely an extremum of the free energy, we may encounter a change in stability or 
a bifurcation point. These transitions and other topological properties of extremal points 
are the subject of Morse theory, which is our first tool. 



4 Namely, we are interested in the topology of the curve of solutions, and their stability, rather than its 
precise location in configuration space. 
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We start by reviewing some features of Morse theory. Consider a function F : M n — > R. 
For each extremal point p where VF = critical groups and a Morse index are defined. 



Denoting F(p) 
homology 



F y = {x : F(x) < y}, the critical groups are given by the local relative 



C*(F,p)=H»(F r ,F r \{p}), 



(3.3) 



where both sets are intersected with a small enough neighborhood of p. The Morse index 
is the (Euler) index of this homology 



^MorseGM = EH 9 dim(C,(F,p)). 



(3.4) 



For the generic case of a non-degenerate Hessian dijF(p), n the number of its negative 
eigenvalues, is called "the Morse type", the critical groups are C n = "Z and all the others 
vanish, and the index is 

J Morse = ( 3 - 5 ) 
For a more general extremal point the index is given by considering the function VF/ \VF\ 
defined on a small sphere that surrounds p as function from S ra_1 to itself and as such it 
has a winding number w and one defines 



Morse 



w{VF(p)). 



(3.6) 



One verifies that the previous definition is indeed a special case of the latter. 

One of the basic results is that if we take an arbitrary sphere S C W 1 with no extremal 
points then 

W (VF(5))=E%orse(p) ( 3 - 7 ) 
p 

where the sum extends over all the extremal points inside S. 

As a parameter is changed transitions may 
occur but the Morse index (and other quanti- 
ties) stay invariant. The most generic tran- 
sition is the annihilation (or creation) of a 
n, n + 1 type pair. For example in Id it is the 
creation of a maximum and minimum pair, 
and can be modelled by the equation F(X) = 
fx X + A 3 (see figure 0) . The more general case 
can be gotten by adding n spectating negative 
modes. 

Out of this primitive vertex one can cre- 
ate more elaborate vertices. The most generic 
vertex for an even F (figure |8|) can be mod- 
elled in Id by F(X) = X 2 — A 4 . It describes 
a solution with n = 1 bifurcating into 2 solu- 
tions with n = 1, and one stable solution (one can think of this vertex as the primitive 
vertex glued to a free line with n = 1). 




Figure 7: The most generic phase transition 
- first order with a cubic term. Stable (unsta- 
ble) phases are denoted by solid (dashed) lines, 
and an arbitrary number of spectating negative 
modes n may be added. 
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Figure 8: The generic phase transition with even free energy, (a) is the phase diagram correspond- 
ing to the free energy F depicted in (b) while (c) is the diagram that corresponds to — F (and — fi). 



It can be seen from these examples that it is useful to think about Morse theory as a 
conservation rule for extremal points under deformations, namely an extremum of type n 
can only be "annihilated" by a vertex with its "anti-particle", an extremum of type n + 1. 

3.3 A suggested phase diagram 

Let us apply the previous discussion to find consistent ways to complete the phase diagram 
in figure (|3|). From the Morse theory "conservation rule" it is now apparent that this 
diagram has a serious problem - if one assumes that a stable non-uniform string phase 
exists for fj, —* 5 there is an inconsistent non-conservation of phases between fj, — > and 
/j, — > oo. The two "completions" of figure (|9|) exemplify the problem: diagram (a) has an 
inconsistent vertex between stable phases, while diagram (b), though consistent, retains 
the paradox of Horowitz-Maeda || intact since the non-uniform phase will have to decay 
now into a BH. 

In order to solve the problem there must be either a new phase at the asymptotic 
regimes fi — > (decompactified 5d) or \i — > oo (4d) or a runaway phase (a phase which 
runs to infinity in field space at finite [i) 6 . Due to these difficulties, and since the stable non- 
uniform black-string was not constructed yet neither analytically nor numerically, I choose 



° which is in spirit of |5j as we discuss shortly. 

6 A11 of these options are problematic, and I speculate here on theorems that could eliminate them. It 
is believed (see also Q section 3.7 and a reference therein) that for fj, — * oo (4d) the only solution is the 
uniform string. One way (suggested to me by J. Maldacena) is to show that as /i — » oo the regime of validity 
of the perturbative analysis around the 4d Schwarzschild solution grows indefinitely. Another way would 



be to compute the index from (3.6). Since the field of metrics is infinite dimensional a winding number 
cannot be defined, but as long as our critical points have finite n it is plausible that the infinite dimension 
could be regularized by a UV cutoff. New localized solutions in decompactified 5d, on the other hand, are 
forbidden by a generalization of the uniqueness theorem EjfJ. 
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Figure 9: Disfavoured completions of the phase diagram. Diagram (a) is inconsistent as it does not 
respect the Morse index at the vertex of the black hole and the non-uniform black string. Diagram 
(b) is consistent but the non-uniform string disappears for small \x and the paradox of Q moves to 
the second transition. 



to eliminate it altogether and suggest the following phase diagram (figure 10), which is the 
simplest option in my opinion. In the suggested phase diagram an unstable non-uniform 
string emanates from the black string at [lql using the even vertex (|8j) and then simply 
annihilates the black hole phase, using a generic vertex ([?]). The latter transition will be 
called a "merger transition" and its fj, will be denoted by //merger- Here the point of phase 
transition is identified with the point of horizon topology change, but a priori these two 
points are independent. Admittedly, other consistent phase diagrams can be imagined, 
but this is the simplest one in my opinion. Therefore the challenge now is to seek some 
confirmation, in which process we will find that changes are due in the neighborhood of 
the merger transition. 



4. Topology - prediction 

Let us use another qualitative tool at our disposal namely the topology of the solutions. 
The black string is contractible to S 2 x S 1 by contracting the (r, t) plane to a point, hence 

string ~ S 2 x S 1 (4.1) 

where the equivalence denotes homotopic (and homological) equivalence. 

The 5d BH is a direct sum of the K 3 x S| x S| background with the 5d Schwarzschild 
solution which is topologically a disk bundle over S 3 . This topology has the following non- 
contractible cycles - an S 1 along the z axis, an S 3 along the horizon, and a less obvious S 2 
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Figure 10: (Left) A suggested phase diagram, which will be refined further in the area of the 
merger transition. \x is the dimensionless parameter, and A is an order parameter (measure of 
non-uniformity) . Dotted lines denote transitions - a first order transition at fix and tachyonic decay 
from the two other points. (Right) Curves of the free energy are given for several /i values, solid 
lines denoting critical graphs for F. 



(see figure 11) which is made out of the segment of the z axis connecting the north and 
south pole, together with the t fiber, which degenerates on the horizon to complete the 
sphere. 

We see that the two solutions have different topologies since the non-contractible 3- 
cycles have a different topology. This is a much stronger statement than the change in 
horizon topology. Usually one thinks about the configuration space of GR as being dis- 
connected between different topologies, namely one is given a manifold on which to solve 
Einstein's equations together with some boundary conditions, but the topology of the man- 
ifold is part of the data, rather than the dynamics. On the other hand we know already 
about several topology change transitions in 6d or higher such as the flop and the conifold 
which occur at finite distance in moduli space. Therefore we continue by assuming that 
such a topology change occurs in our 5d system as well, and seek a local model for it. 



Analyzing the topology near the merger transition (see figure 11) we see that the 
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topology of a 4d boundary away from the wall (which hence remains constant during the 
transition) is S 2 , ± x S 2 zt - the line in the (r, z) plane combines with the t circle fiber to 
give S 2 z t over which the S 2 , , is fibered at each point and the fibration can be seen to 
be trivial due to the decoupled form of the metric 2.2. In the black-hole phase the S 2 , , 



shrinks completely, while in the string phase the S 2 zt shrinks. So the topology change 
can be modelled by the "pyramid" familiar from the conifold transition, and the singular 
solution is the cone over S 2 x S 2 . 





Figure 11: The merger transition. Shaded regions are inside the horizon and the dashed line is 
the base. The singular configuration is a cone over S 2 x S 2 . 

Topologically the transition can be modelled by the equation 

\x\ 2 - e = \y\ 2 (4.2) 

where f S I 3 is a 3-vector and , y = (r,t,z) combines the z coordinate nicely with (r, t), 
the origin of z being at the merger point, and the origin of (r, t) on the horizon. For e > 
we have the string, while for e < it is the black hole. We view the realization of the 
transition within the class of Ricci-flat metrics as a challenge for our picture. 

5. A local topology change - confirmation 

It is required to find Ricci-flat metrics with 50(3) x U(l) isometry which realize the 
topology change through the cone over S 2 x S 2 . 

5.1 The cone 

We start by looking for a Ricci-flat metric for the cone itself, ds 2 = dp 2 + p 2 rfs| 2xS2 , 
which is equivalent to finding an Einstein metric on S 2 x S 2 with a positive Einstein 
constant 3 (the unit S 2 has Einstein constant 1). I find that the most general such metric 
is 3ds| 2xS 2 = d\ 2 + e 2c sin 2 (x) da 2 + d6 2 + sin 2 (6>) d(j) 2 , where c is an arbitrary constant 
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which will be set to zero to avoid a unwanted singularity which persists away from the 
apex. So the geometry is just a product of two round S 2 's with the angle of the cone 
chosen to assure Ricci flatness 

ds 2 = dp 2 + ^-(dn 2 ^ + dnl a ). (5.1) 

Note that the isometry of the local solution enhances to SO (3) x 50(3). 

Let us verify that this is the only possibility. The 50(3) x U(l) isometry generates 
homogeneous 3d orbits [50(3) x U(\)]/U(l) which are classified by an integer k giving 
rise to S 3 /^^ - a circle bundle over S 2 with monopole charge k (with the Wjj, acting 
along the Hopf fibration). Then S 3 /^^ must be fibered over an interval (or circle but 
that would give a trivial fibration and the wrong topology) and "sealed off" at the end- 
points to make a compact manifold of topology S 2 x S 2 . The possible metrics are ds 2 = 
dx 2 + e 2 A {dO 2 + sin 2 (0) d(p 2 ) + e 2B (da -k cos(9)d<p) 2 . Moreover the discrete TL 2 isometry of 
our solution for reflection of the z axis, translates into the functions A = A(x), B = B(\) 
being even. If the function e 2B has zeroes at the end of the interval then we get the 
Hirzebruch spaces, - fibrations of CP 1 over CP 1 . Therefore the only way to get the 
product topology S 2 x S 2 is to start with k = 0. Then the Einstein equations determine 
e 2A = ci, e 2B = C2 sin 2 (x), as claimed in the previous paragraph. 

5.2 The deformed cone 

We now proceed to look for Ricci-flat solutions over the deformed cone. A priori I do not 
know of any obstruction to do that: in 2d the deficit angle of the cone would impose some 
non-zero curvature (and Ricci) due to Gauss-Bonnet, but that does not apply here. 
The ansatz is 

ds 2 = dp 2 + e 2a (d9 2 + sin 2 (6>) d(f> 2 ) + e 2b (d X 2 + sm 2 { X ) da 2 ) (5.2) 

where a = a(p), b = b(p) and we retained the enhanced SU(2) x SU(2) isometries. 
The equations are 

a" = e - 2a -2a' 2 -2a'b' 
b " = e -2fe _2b' 2 -2a' b' 

= a' 2 + a" + b' 2 + b" (5.3) 

One verifies that the 3 equations for 2 unknown functions imply a first order constraint 

= a ,2 + 4a'b' + b ,2 -e~ 2a -e- 2b (5.4) 

which is consistent with the second order equations, and reduces the number of required 
initial conditions from 4 to 3. 

The cone solution corresponds to ao = &o = log(p) — log(3)/2. Let us linearize first the 
equations around this solution. Denoting 

a = ao + a b = b$ + (3 (5-5) 
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we find the linearized equations 



6 (4 + 2) , 2 
a = — t a a (3 

P l P P 

0I= 6 (4 + 2) 2_ a , 

P Z P P 

= a" + -a' + 13" + V (5.6) 
P P 



Due to the symmetry in exchanging a and (3 it is easy to diagonalize by choosing a 
a + (3, a~ = a — f3. The solutions are 



„+ 



a+ = c+p s i +4/2 = -1, -6 ^ ^ 

a~ = c7/p s i +c^p s 2 s ~ 2 = (— 3±i\/l5)/2 

where c^ 2 are arbitrary constants. The interesting modes are the modes of a~ , since the 
solution with s + = —6 is forbidden by the first order constraint, while the s + = — 1 is simply 
a translation in p (a diffeo) which had to be there since the equations are p independent. 

We now proceed to solve the non- linear set of equations (|5.3| ). Clearly the linear 
approximation is expected to hold only for large p when the corrections are small. We are 
looking for a solution where for p — > e 2a will combine with dp 2 to make a smooth 3d 
space, while e 2b approaches a finite limit. Therefore we expand around the singular point 
p = taking the following ansatz 

a(p) = \og{p) + ai p 2 + a 2 p 4 + a 3 f + a 4 p 8 . . . 

b(p) = log(po) + h p 2 + b 2 f + b 3 p 6 + +6 4 P 8 ■ ■ ■ (5.8) 

where p := p/po, and we hope the ansatz is self consistent and allows to solve for the series 
coefficients a^, hi order by order. Indeed, we find using Mathematica [^] the equations to 
be solvable and 



ai = -1/18 --5.6-10- 2 6i = l/6 -1.7-10- 1 

a 2 = 67/8100 ~ 8.3 • 10~ 3 b 2 = -13/540 ~ -2.4 • 10~ 2 

a 3 = -4031/2679075 ~ -1.5 • 10" 3 b 3 = 191/42525 ~ 4.5 • 1(T 3 

a 4 = 163777/535815000 ~ 3.1 • 10" 4 64 = -33391/35721000 ~ -9.3 • 10~ 4 



(5.9) 



After finding these solutions I was disappointed to find that their existence was al- 
ready discussed in the mathematics literature p7| , p8[ where similar Ricci-flat metrics of 
cohomogeneity one are described which approach a cone as a limit, and that the initial 
value problem was discussed in [^{J, although these papers do not consider going through 
the cone to achieve a topology change. 

5.3 The cone and the action 

Let us estimate the dependence of the action on the smoothing parameter pq. Since the 
metrics are solutions for all p the bulk term vanishes and only the boundary term con- 
tributes to the action. We saw that at large distance distances the corrections to the 
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metric behave as h ~ (p/ Po)~ 3 ^ 2 ■ Thus I ~ /i 2 ~ p^/G^. The proportionality factor could 
include a log(pi/po) factor as well as an arbitrary multiplicative constant, where p\ is a 
long distance cut-off on the validity of the local model, p\ being of order L which is the only 
length scale at /^merger- The multiplicative constant can and in general will be different 
on the two sides of the transition, producing a high derivative "kink" in /, namely a phase 
transition point. 

Despite the success of describing the a Ricci-flat topology change using the cone over 
S 2 x S 2 it turns out that this cone is unstable? 

Let us first verify this result and then discuss the implications for the phase diagram. 
The ansatz for the perturbation is 

n 2 

ds 2 = dp 2 + P — {e 2e dn 2 d!<j} + eT 2e dQ 2 Xia ) . (5.10) 
o 

Namely, the unstable mode, e(p), is a change in relative size of the two spheres. The action 
for this mode is 

1 = I (t ) 2 pi [ 12(cosh(e) - 1)1 p2 - e ' 2 ] ■ (5 - n) 

Since F = —1/(3 we see that e{p) indeed has a negative mass-squared (actually, due to the 
volume factor it is not evident yet that the mode is tachyonic as seen below is detail). 

The instability is alarming at first, but it can merge nicely in the picture. The point is 
that if the cone is cut-off both at small and large p then it becomes stable. The cut-off at 
small p is supplied by po while at large p the local cone description ceases at a length-scale 
pi of the same order as L or the Schwarzschild radius. 

Let us estimate when should this instability appear. The quadratic part of the action, 
throwing out overall constants is 



/ 



I ~ p dp 



— e 2 -e' 2 

i 

P 



(5.12) 



Looking for a marginally stable mode (a zero-mode) we should solve the equations of motion 
e" + 4e' / p + 6e/ p 2 = 0, whose solutions are e{p) = Yli=i c i P Si where the exponents are 
given by s\^ = (—3 ± i\/TE)/2, which are naturally the same exponents that appeared in 
the linearized analysis (eq |5.7| ). Hence if Dirichlet boundary conditions are put at po,p± 
marginal stability will occur at 

pi/po = exp(vr/Im(s)) ~ 5.1 (5.13) 
Alternatively, the action ( 5.12j ) can be recast in canonical form 



I ~ J dp 

through the change of variables p -1 = 3 p s . 



■ e z - e z 



3p 2 



(5.14) 



7 I thank I. Klebanov for suggesting this, and M. Rangamani and E. Witten for sharing this unpublished 
result Boll. 
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6. Generalization and critical dimension 



Let us consider generalizations of the previous section, section ||. Instead of starting with 
a background K 4 x S 1 we could have taken IR^ -1 x S 1 with d > 5 creating a transition 
in which the S d ~ 3 of the string horizon shrinks, and is replaced by an S 2 connecting the 
north and south pole of the BH along the z axis. This would lead us to study cones over 
S 2 x S d - 3 . 

Here we choose to study more generally cones over S m x S n because we can do this 
more general case for the same price. It would be interesting to find which of these and 
of other general cones can be produced by replacing S 1 by a more general compactifying 
manifold. 

It will turn out that all the important properties of these cones will depend only on 
the total dimension d = m + n + 1, and so it is natural to speculate that this property is 
more general. 

The Ricci flat cone is given by 

o 2 

ds 2 = dp 2 + ((m - 1) dn% m + (n - 1) dJ2|») , (6.1) 

where the constants are required in order to get the correct Einstein constant for the base, 
namely d — 2 (S m has Einstein constant m — 1). 
The ansatz for a deformed cone is 

ds 2 = dp 2 + e 2 a(p) dn 2 S m + e 2 b{p) dn 2 S n (6.2) 

and Einstein's equations are 

a" = ( m - l) e - 2a - ma' 2 - na' b' 
b" = ( n -l) e - 2b -nb' 2 -ma'b' 

= m {a' 2 + a") + n {b' 2 + b") . (6.3) 

One verifies again that the equations imply a first order constraint which is consistent with 
the second order equations, and reduces the number of required initial conditions to 3. 

The cone solution corresponds to a = log(p) + log(m — l)/2 — log(d — 2)/2, bo = 
log(p) +log(n — l)/2 — log(d— 2)/2. Let us linearize first the equations around this solution. 
Denoting 

a = ao + a b = bo + (3 (6-4) 
we find the linearized equations 

2 (d — 2) 2 m + n . n . 
a = ^ a a p 

P P P 

0fl = _2(d-2)p_2n + m rn a> 

P 2 P P 

= m (a" + -a') +n((3" + (6.5) 
P P 
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The equations are diagonalized by choosing a 1 
are 



a 

of 



cl p i + cTJ p 
Ci p Sl +C2 9 s 



'1,2 



'1,2 



ma + n (3, a 
-1, -2(d- 2) 



^(-l±h/l-^2 



a — (3. The solutions 
(6.6) 



where cf 2 are arbitrary constants. Again the interesting modes are the modes of a~, 
since the solution with s + = — 2(d — 2) is forbidden by the first order constraint, while 
the s + = — 1 is simply a translation in p which had to be there since the equations are p 
independent. 

Note that d = 10 is a critical dimension for the discriminant. For d —* 00 the two roots 
are real and approach —{d — 2) and —4. 

Looking for a smoothed cone we repeat the expansion of a(p), b(p) near p = as in eq 
flO|), and find 

— n (n— 1) 



ai 6m(m+l)' 

(n— l) 2 n (l8m 2 — 3n+17mn) 

° 2 ~~ 180 m 2 (m+1) 2 (m+3) 

(n-1) 3 n (405 m 4 +30 n 2 -119 m n 2 +27m 3 (15+32 n)+m 2 n (270+461 n)) 

a 3 - 



bi 
b 2 



("-!) 
2 (m+1) 

-(n-1) 2 (3+3m+2n) 
12 (m+1) 2 (m+3) 



5670 m 3 (771+1)^(771+3) (m+5) 
90m(m+l) T "(15+8m+m 2 ) 



^ _ (n-1) 3 (l5m :i +n 2 +m 2 (45+24 n)+m (30+30 n+11 n 2 ) ) 



(6.7) 



Let us turn now to the stability analysis. The fluctuations we want to analyze are 

„2 

1 1 ... i ) . 2<n,> ../o- i / ., ii,--' » ,?r £ 



ds 2 = dp 2 + ( (m - 1) e 2 £/m d^l™ + (n - 1) e" 2 e/n dOln ) . (6.8) 



d-2 

The action for this mode is 



cp 



4 fm (d - 2) e~ 2e / m + n (d - 2) e 2e / n - (d - 1) (d - 2)) - (- + — ) e' 



(6.9) 

where the constant is given by c = O m O n (m — l) m / 2 (re — l) n / 2 (d — 2)~( d-1 )/ 2 and S7 m = 
(to + 1) 7r( m+1 )/ 2 T((m + l)/2 + 1) is the volume of S m . The quadratic part of the action, 
disregarding overall constants is 



„<2-l 



dp 



2 (d - 2) g2 £/ 2 



(6.10) 



The zero mode equation is e" + (d — 1) e'/p + 2 (d — 2) e/p 2 = 0. The solutions are e(p) 
Z)i=i c i P Si where the exponents are given by 



S\,2 



( _1 ^ l \l 77=2 



1) 



(6.11) 



Again these are the same exponents we got in the linearized analysis. We see that 
the field e(p) is unstable only for d < 10 while for d = 10 it is marginal in the quadratic 
approximation. Let us tabulate 



d 


5 


6 


7 


8 9 


10 


Im(s) 2 


15/4 


4 


15/4 


3 7/4 





exp(7r/Im(s)) 


5.06 


4.81 


5.06 


6.13 10.7 


oo 



(6.12) 
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In general 



4Im(s) 2 = 16- {d-§f 



(|34|], a closely related formula appeared in |35| eq. (34)) and 



(see figure p"^ ) 



Re(s) 



-2Re [s! 



-(d-2)/2 



for d < 10 



-(d-2)/2(l± A /l- g^) ford>10 



(6.13) 



(6.14) 




Figure 12: A plot of — 2Re(s_) as a function of the dimension d. 
Alternatively, the action ( |6.10| ) can be recast in canonical form 

'~/*[oni?''"''1 <fU5) 

through the change of variables p^ 1 = (d — 2) p d ~ 2 . 

7. Summary and Discussion 

7.1 Summary of suggested phase diagram 

Let us recapitulate the suggested picture for the transition which emerges from our discus- 
sion (figure |l0|). The transition has a typical behaviour for a first order transition, in the 
presence of an even free energy, the only stable phases are the black-hole and the uniform 
black-string and there is an unstable non-uniform black-string phase. Starting with the 
string and reducing p (string — ► BH) one reaches first fi\ the point of first order transition 
to the BH phase. However, since the transition occurs via tunnelling due to statistical 
and/or quantum fluctuations the string is still meta-stable. If we continue to lower \i we 
reach the Gregory-Laflamme point, pgl, where a marginally tachyonic mode appears and 
the system decays immediately. For lower p the string is still a solution, though unstable, 
and we encounter "copies" of the transition at \i — ► /j,/k, k £ N. 
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Going in the opposite direction, BH — * string, one starts with a 5d black-hole at small 
\i. As jj, is increased one encounters the first order point fj,\ beyond which the BH is only 
meta-stable and later one reaches a perturbative instability. Altogether the system traces 
a "hysteresis" curve typical of a first order transition. 

In our first suggestion the perturbative instability was thought to occur at /x me rger 
which is where the black hole phase terminates. However, it turns out that the cone 
over S 2 x S 2 is unstable (@ and related results in for AdS a x M p x M 7 _ p , g for 
AdSj x S 2 x S 2 and |33| for AdS 5 x S 3 x S 2 ), and a dramatic dependence on the dimension 
is found in which cones over products of spheres have an instability for d < 10 which 
disappears for d > 10. Hence for d > 10 I view the original suggestion (figure [R]) to be 
viable, and a small feature must be added, namely a high derivative kink (2rd derivative 
kink in 5d) to the curve /x = //(A) at the merger point (see subsection For 5 < d < 10 
on the other hand the picture must be more interesting and we must amend the phase 
diagram at the vicinity of the merger transition expecting the perturbative instability to 
occur before merger, roughly when ri/ro — D ~ 5 where r\ represents the radius of the 
spheres at the large length cutoff of the cone - at the outer edge of the cone, and r is the 
radius of the minimal S 2 , a. 

It is tempting to speculate further about the form of the phase diagram in the vicinity 
of the merger transition. Consider decreasing r$ while r\ stays roughly constant. The 
resolved cone continues to be a solution as we saw, but each time we cross a threshold 
r i/ r o — D k , k 6 N a new marginally tachyonic mode emerges. The situation is just like 
moving along the string phase - each time a critical point is reached two phases emerge 
from this point. A detailed local analysis is required in order to determine whether these 
new phases are stable or not (the analog of the analysis of § at the GL point). One 
possibility is that the phase diagram near the merger point will look like repetitious copies 
- see figure 13(a). Since the cone is self-similar and the new phases appear again and again 
in a self similar way it is tempting to speculate about a more involved possibility with a 
fractal phase diagram and fractal physics near the merger point (see figure |l3|(b)), in which 
the constant D will be related to the fractal dimension. This will be in tune with the 
numerical discovery of scaling behavior in General Relativity |p9| . 

The topology changing transition does not occur for d < 5. Physically it is because 
that would require analyzing black holes in 3d and lower which do not have regular horizons 
but rather induce deficit angles, while mathematically it is because both sphere factors in 
S m x S n must have positive Einstein constant and so m, n > 2. For 5 < d < 10 the picture 
stays similar to d = 5. 

However, for d > 10 the previously tachyonic function stabilizes due to a competition 
in the action between the volume factor and a certain mass-like term. Closely related 
dependence on the dimension was observed in geometric studies of minimal cones p6|; 



in construction of cohomegeneity one Ricci flat metrics (27], 28]; in the stability of AdS 



compactifications over a product space |5g] and J37[ (this result could be related to this 



one by the cone- AdS relation |38| , though here D brane probes do not exist as such) ; in [ 34 ] 
and in the stability of Schwarzschild black holes in a asymptotically conical background 
1 35 ] . However, here they arise in a natural dynamical context for the first time. I can 
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Figure 13: Possible behaviors around the merger point for 5 < d < 10. (a) minimal change, (b) A 
fractal diagram. The thin dashed line separated the black hole phase on the right from the unstable 
string to the left. The number of negative modes for each phase is suppressed here. 



see no connection between the lOd criticality observed here and the critical dimension of 
superstring theory, and so if such a connection existed it would be very interesting. 

A phenomenon that could remind one of the Gross- Witten transition as claimed in Q 
was found, namely, that the free energy probably has a kink on the wall in configuration 
space where a topology change takes place (subsection |5.3| ). In 5d the jump is in the 
3rd derivative (or higher), and in other dimensions it is at least — 2Re(s_) (figure 12). 
Moreover, a high order phase transition from black hole to non-uniform black string seems 
to occur at the merger point, but this is not a real phase transition, unless the emanating 
non-uniform black string is stable (which I do not expect). This issue should be revisited 
when the phase diagram in the vicinity of merger transition is better understood. 



7.2 No Cusps 

The cusps on the horizons which were discovered in numerical simulations of black hole 
collisions (see and references therein) were not seen in the local model which has smooth 
minimal spheres. That seems to indicate that the cusps are velocity dependent, vanishing 
for a static configuration. It would be interesting to test this by measuring numerically the 
cusp for various relative velocities keeping equal separations. Another possibility is that 
the local model is not sensitive enough to detect the cusps. 



7.3 Topology change 

The term "topology change" in General Relativity is being used with several different 
meanings, for instance the topology of the whole space-time could change as a parameter 
is varied or the topology of a spatial slice could change as time progresses. The topology 



change here will be of the first kind, evolving with a parameter, just like the flop [40] 



and conifold [41] transitions for 6d Calabi-Yau surfaces. Expected implications for time 



evolution are briefly discussed below, and we should note that classical no-go theorems 



|4q , 46] are not valid due to the presence of the singularity (somewhat in the flavor of the 



possible classical topology changes of 



- 22 - 



Let us compare the merger transition discussed here with the conifold transition [41] 



and the "Gr2-conifold" 44]. The merger transition is in m + n + 1 dimensions and 
centers around the cone over S m x S n (minimum 5d) while the conifold (GVconifold) is in 
6d (7d) involving the cone over S 3 x S 2 (S 3 x S 3 ). Despite the similar topologies and Ricci 
flatness of the geometries they differ in supersymmetry and metric: whereas the previous 
topology changes had Killing spinors the merger cone does not, and indeed it is not even 
stable. Accordingly the metric on the spheres differs - the conifold (GVconifold) metric 
is given by the homogeneous space T M = (5*7(2) x SU{2))/U{1) (SU "(2) 3 / 'SU '(2) diag ) 
whereas the merger cone involves the round spheres and in 6d can be written as T±o = 
SU{2)/U{1) x SU(2). 

7.4 Future directions 

I would like to point out the following topics for future work 

• Phase diagram near the merger point. 

• Lorentzian signature 

• Explosive transition 

• Uniqueness theorems (no hair) 

• a' corrections 

• Merger and quantum gravity 

which will be discussed below one by one. 

The suggested phase diagram (figure ^) was found to require corrections near the 
merger point at least for 5 < d < 10. Some possibilities were discussed above in subsection 



7.1 



So far we worked only in the Euclidean signature, and it is an important question to 
understand the Lorentzian. The static solutions can be readily turned to static Lorentzian 
solutions outside the horizon, and one expects an extension into the inside of the horizon 
to exist. However, the singular cone solution at the merger transition will translate into 
a naked singularity on the horizon, at which point the structure of the light cone may 
be different and a problem with the extension of the solution could arise. Even more 
interesting is the time dependent behavior, when the tachyon decays. For instance, one 
could start with an unstable string fi < fiQL with initial conditions that as t — > — oo, 
both the tachyon its time derivative approach zero T, T — ► 0. Since the problem is time 
dependent the apparent horizon does not coincide anymore with the event horizon, and it 
would be interesting to understand the Penrose diagram. 

The hysteresis curve contains two sections of tachyon decay, one starting after the 
Gregory-Laflamme point, and the other close to the merger transition. Since the energy 
emitted during these decays scales as the mass of the black object, it must equal a numerical 
constant times the mass (there are no dimensionless parameters at the transition points). 
Assuming large extra dimensions one finds large fusion/fission explosions [j49f| . 
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The black hole uniqueness theorems ( "no hair" ) , which characterizes a time-independent 
4d black hole by its mass and angular momentum (and possibly gauge charges), does not 
hold as such for d > 4. The 5d rotating black ring solution of |Q is a beautiful counter- 
example where for some values of the mass and angular momentum 3 solutions exist - 2 
black rings and 1 black hole. The black-hole black-string transition provides another ex- 
ample for non-uniqueness above 5d, since at some values of [i there are at least 3 solutions 
- a (perturbatively) stable string, an unstable non-uniform string and the black hole. This 
example is static, but the asymptotic space-time is not Euclidean (though it is flat). It 
would be a shame to loose such a nice principle at higher dimension, and it would be better 
to try and "fix" it. Following this work a natural speculation arises - that the topology of 
the horizon must be specified, and that unstable solutions should not count [50]. 

It would be interesting to find string theoretic a 1 corrections to the cone geometry. On 
the one hand large corrections could appear close to the tip where curvatures approach the 
string scale, but on the other hand the cone topology seems to be essential for the topology 
change, which may protect the cone. Perhaps non geometric modes become important like 
in the conifold transition 

Finally it is natural to wonder about the quantum gravity interpretation of a black 
hole merger, as mentioned in the introduction, and it is hoped that this work makes a step 
in that direction. 



Acknowledgements 

It is a pleasure to thank Lenny Susskind for introducing me to this problem some years ago; 
Gary Horowitz for rekindling my interest and for discussions; John Bahcall for encourage- 
ment; Tsvi Piran and Evgeny Sorkin for collaboration on a related project and comments 
on the paper; Sunny Itzhaki and Juan Maldacena for some enjoyable discussions; Steve 
Gubser, Igor Klebanov, Mukund Rangamani, Edward Witten for correspondence and im- 
portant comments; David Berenstein, Sergey Cherkis, Eric Gimon and other members of 
the IAS group for discussions; and finally the University of California at Santa Barbara 
and Stanford university for hospitality while this work was in gestation. 

The first version of this paper appeared while I was at the School of Natural Sciences, 
Institute for Advanced Study, Princeton and the work was supported by DOE under grant 
no. DE-FG02-90ER40542, and by a Raymond and Beverly Sackler Fellowship. Currently 
my work is supported in part by The Israel Science Foundation (grant no 228/02) and by 
the Binational Science Foundation BSF-2002160. 



A. Comments on a local analysis around the GL critical point 

In this appendix I would like to describe an approach to the determination of the local 
properties of the system at the GL point which is somewhat different then the one used 
in H and is actually standard in the thermodynamics of phase transitions. It requires the 
computation of solutions up to the second order in perturbation rather than third, but the 
required hard computations were not pursued to the end. 
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Start by considering the action functional, namely working in the canonical ensemble. 
At the GL point the quadratic potential has one zero mode, and the objective is to go 
beyond the quadratic approximation, and determine whether this mode is negative or not, 
which are the only qualitative possibilities. This will determine the current goal, namely, 
whether the emerging phase is stable or not, and whether its \i starts climbing down or up. 

The quadratic approximation to the action can be written as 

-h\x,yi\fi] = k £i\x\ 2 + ^7?i||yi| 2 (A.l) 

i 

where x is the marginally stable mode, yt are all the other modes diagonalized such that 
their mass-squared are m 2 , jX = fj, — [iql and fco is a positive constant. One would like to 
trace the emerging solutions x , y i , fJj clS cl function of a perturbation parameter e. Clearly 
at first order only x can be turned on and the first task is to find it explicitly by solving 
the equations of motion for I2. 

Due to the U(l) z isometry of the uniform black-string all modes with z-dependence 
are complex with a specific charge qt such that under a rotation of the z axis by an angle 
9 they transform by 

y l ^e^ e y l (A.2) 

By the construction of x we have q x = 1. Let us use the fact that the action is real and 
invariant under this action. The first consequence is that cubic terms involving only x, x 
(which would imply a first order transition) are forbidden. Therefore we must go to quartic 
order in e. For this we must expand the action to include 

-I 4 = gi yi x 2 + c.c. + hj yj \x\ 2 + k |x| 4 (A. 3) 

» 3 

where the sum is over modes of charge ±2 or which couple to quadratics in x, x and gi, hj 
are the respective couplings. 
One takes 

x = x m e + x m e 2 

Vi = yf ] e 2 (A.4) 

And finds that — (I2 + h) becomes a quadratic form in e 2 , and its signature is determined 
by the sign of 

*-&W-Efrj (A.5) 

i 3 

In practice it is more convenient to find the solution for x, y% up to second order in 
e and check the sign of the resulting action — (I2 + -^4) - a positive mode means that it 
is a minimum and hence a second order phase transition will form, whereas a negative 
mode means that there is already another minimum lower than the present one, and so the 
system goes through a first order transition (see figure ||) . 



If one computes the constant ko in I2 flA.lD one can solve for e = e(/t) and extract 
quantitative information about the initial rates of change of /, x, yi along the new phase 
as a function of ji. 
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